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Abstract We study consistently the pion's static observables and the elastic and 7*7 — > ir° transition 
form factors within a light-front model. Consistency requires that all calculations arc performed within 
a given model with the same and single adjusted length or mass-scale parameter of the associated 
pion bound-state wave function. Our results agree well with all extent data including recent Belle 
data on the 7*7 — > n° form factor at large q 2 , yet the BaBar data on this transition form factor 
resists a sensible comparison. We relax the initial constraint on the bound-state wave function and 
show the BaBar data can partially be accommodated. This, however, comes at the cost of a hard 
elastic form factor not in agreement with experiment. Moreover, the pion charge radius is about 40% 
smaller than its experimentally determined value. It is argued that a decreasing charge radius produces 
an ever harder form factor with a bound-state amplitude difficultly reconcilable with soft QCD. We 
also discuss why vector dominance type models for the photon-quark vertex, based on analyticity and 
crossing symmetry, are unlikely to reproduce the litigious transition form factor data. 

Keywords Neutral Pion • Axial anomaly • Form factors • Light-Front Field Theory 



1 Introduction 

The transition 7*7 — > ir° has attracted considerable attention with the experimental findings of the 
BaBar Collaboration [1]. Indeed, while this new data set appears to agree with earlier experiments 
on a domain of squared-momentum transfer below Q 2 = —q 2 < 10 GeV 2 [2; 3], the data points at 
larger Q 2 values remarkably exceed the prediction of perturbative QCD (pQCD) in the asymptotic 
limit [4; 5]. On the other hand, a most recent measurement by Belle [6] appears to corroborate the 
pQCD limit. The process 7*7 — > tt° is by itself of highly phenomenological and theoretical interest. 
If the entire domain of experimentally explored squared-momentum transfer is to be described within 
a unique theoretical framework, it must account for the mainly nonperturbative phenomenon of the 
Abelian anomaly and functional behavior of perturbative QCD. This is the challenge for any model 
employed in a self-consistent way. 

Herein we present the results of a consistent treatment of the pion's static features as well as of 
the electromagnetic and 7*7 — > tt° transition form factors within one particular model, namely the 
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light-front model introduced in Ref. [7; 8] and refined in later studies [9; 10; 11; 12; 13; 14; 15]. This 
work differs from several recent attempts [16; 17; 18; 19; 20; 21; 22; 23; 24; 25; 26; 27] to reproduce the 
large-Q 2 BaBar data in spirit and in aim: we do not attempt to build a phenomenological model that 
yields a good fit to the BaBar data. Rather, within the framework of the present model, we insist on 
a consistent and simultaneous treatment of all extent data, not merely the 7*7 — > ir° transition form 
factor. 

We first calculate the static properties, the pion decay constant and electric charge radius. The pion 
decay constant, f„, serves to adjust the unique parameter that enters the bound-state wave function, 
r nr , and introduces a length (or mass) scale. This immediately fixes the root- mean-square (rms) charge 
radius of the pion whose known experimental value, (r 2 ) 1 ' 2 = 0.672 ± 0.008 fm [28], can be compared 
with. Next, with this parametrized bound state function, we calculate the leading contribution to 
the electromagnetic form factor, F^IQ 2 ), and to the transition form factor, F 17r (Q 2 ), in the impulse 
approximation on the light cone. The main freedom in these calculations is, of course, the choice of 
bound-state wave function or in other words, the associated Bcthe-Salpetcr amplitude. 

We devise a wave function which, in interplay with our light-front quark model, reproduces the 
experimental space-like tail of F 7r (Q 2 ) and thus the asymptotic limit of the product Q 2 F V {Q 2 ) — > 
167r/ 2 a s (Q 2 ) for Q 2 — s- 00 [4; 5]. It should be stressed that the decay constant and charge radius are 
insensitive to the asymptotic form of the wave function. This is the reason for the success of quark 
models or related contact-interaction models in the calculation of the pion's static properties below 
a mass scale M 2 [29]. However, probing the pion with Q 2 > M 2 leads to marked deviations from 
experiment when the contact-interaction model is treated self-consistently with proper regularization 
[30; 31; 32]. In particular, due to the lack of a running mass function, models based on constituent 
quarks must necessarily introduce a phenomenological bound-state amplitude that mimics both soft 
QCD effects and the asymptotic behavior of the associated pion distribution function, (p v (x). Therefore, 
this process also allows us to scrutinize and improve the description of the pion bound-state function 
in relativistic light-cone models. 

To summarize, we first describe the light-front formulation of the 7*7 — > ir° and introduce the 
parametrization of the pion bound-state wave function in Sec. 2. The values we obtain from this 
paramctrization for F J7r (Q 2 ) and F 7T (Q 2 ) as well as the pion charge radius and decay constant are 
discussed in Sec. 3. We resume in Sec. 4 where we relax the condition on the length scale parameter 
and observe the implication on the functional behavior of -F 77r (Q 2 ), in particular whether the asymptotic 
prediction of perturbative QCD is still verified. We also check the values obtained in this case for (r 2 ) 1 / 2 
and f„ and modifications occurring in F W (Q 1 ). In Sec. 5 we wrap up with a qualitative discussion of 
the 7*7 — > 7T° form factor and its possible evolution upon the inclusion of vector resonances in the 
deep space-like region. 



2 Pion form factors on the light-front 

2.1 Preliminaries 

The geometry employed, for example by the BaBar Collaboration, to produce the ir° consists of an 
untagged almost real (fcf ~ 0) photon scattered at small angle from the collision axis and a tagged 
electron emitting a highly off-shell photon with space-like momentum k\ . The two photons eventually 
create the neutral pion via the Abelian anomaly. 

The relativistic approach to the wave function based on constituents quarks is possible due to the 
absence of pair-creation processes on the light front. This salient feature arises from the particular 
choice of the light- front coordinates [33; 34; 35], defined by x + = x° + x 3 , in which the center-of-mass 
is readily separated [36]. In essence, the impulse approximation depicted in Fig. 1 involves a bound 
state and an internal loop momentum. The integration over the loop momentum is first performed in 
the convergent light-front energy, k~ , after which the pion's wave function emerges naturally within the 
three-dimensional remainder of the integral and depends on the + and _L component of the qq relative 
momentum [7; 8; 9; 10; 37]. We stress that in defining a model for this light-front wave function, the 
convergence of the fc~ integration is crucial: it is the case for the good component of the current [33; 34] 
and also for the diagram that describes the weak decay constant. 
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Fig. 1 The leading diagram corresponding to the impulse approximation of the process 7*7 — > tv . The lines 
with arrows denote constituent light quarks; the solid blobs represent bare photon-quark vertices and the filled 
oval depicts the qq bound-state wave function of the pion. 



The aim of this paper is to study models of light-front wave functions in computations of the 
electromagnetic and transition form factors of the neutral pion with simultaneous verification of the 
related weak decay constant and charge radius. We make use of two distinct model wave function: 

i) a Gaussian model [38; 39]; 

ii) a hydrogen- type wave function [40]. 

Both wave functions feature a length or mass scale parameter, r nr , whose value is of the order of 
^Iqcd- We fix this scale with a fit to the pion decay constant, all other observables follow from this 
parametrization without exception. We compare the results for the form factors obtained with these 
two wave functions in Sec. 3. 



2.2 Formulation of light- front amplitudes and form factors 

In the following, we formulate the 7*7 — > ir° transition form factor in the light-front approach. 
The coupling of the pscudoscalar qq pair to the pion is expressed through the following effective 
Lagrangian [7; 8], 

C^=-z^-q^rq, (1) 

where M is the constituent-quark mass, /„■ = 92.4 MeV is the weak decay constant of the pion, it and 
q are respectively the pion field and quark wave functions. The units are chosen such that H = c = 1. 
This coupling can be thought of as the leading term of the full pseudoscalar Bethe-Salpeter amplitude. 

The amplitude for the electromagnetic process 7*7 — > tt° in the impulse approximation is diagram- 
matically represented in Fig. 1 and given by the tensor, T^ v ', which contains two contributions due to 
bosonic symmetrization of the amplitude: 

T^(h,k 2 ) = t liv (k 1 ,k 2 ) + t pv (k 2 ,ki) . (2) 

The tensor t fiu (kx, k 2 ) is obtained after taking the traces in spinor and flavor space [41]: 

4 M 2 p 

where I(kf) is the scalar loop integral, 

7 ^ = / (2^0* ((fc 2 - fc) 2 - M 2 + k) (fc 2 - M 2 + i€){{K - kf - A/ 2 + ie) ' ^ 

Here, N c is the number of colors, k^ = k\ + k 2 is the ir° momentum, k\ = q is the space-like momentum 
transfer, k 2 = is the on-shell photon and eo and M are the unit charge and constituent quark mass, 
respectively. The factor 1/3 in Eq. (3) stems from the flavor trace. 



4 



After transformation to light- front variables, k^, k + = k° + k 3 and fc~ = fc° — k 3 , we first integrate 
over k~ , the light- front energy. The reference frame is chosen such that q + = q~ = and momentum 
transfer is transversal, q±, which is always possible due to the space-like character of q. In this coordi- 
nate system, and after factorizing k + . k+ — k + and k^ — k + in the denominator, the integral in Eq. (4) 
becomes 

1 



_ , (k a -k)j+M'- t6 \ /,_ _ ,_ (k.-k)» +M'- t6 
K 2 K K7r K TIT 



(fcj-fc+) / V (fei-fe+) 

where fc+ = A:^ and fc~ = k^ in this reference frame. 

The integration in k~ , performed using Cauchy's theorem, is convergent. The position of the poles 
in the k~ complex plane depends on the value of k + . Since fc+ = k%, this integral is non- vanishing in 
only one region of k + . The set k + < and k + > k+ does not contribute to the integral because the 
three poles in there have imaginary parts with same signs. The only contribution comes from the region 
where the + component of the quark momentum cannot exceed the pion momentum: < k + < fc+. 
This corresponds to an on-shell quark. The result is 

r, 2n = - % f dx d2k ± 1 / fi N 

^ ' 2(2tt) 3 J x(l - x) 2 ( k - k + _ fcj+jg _ (k 3 -k)»+Af' y fc _ fc+ _ kl+AP- _ (k T -k)j + M2 \ ' ^ > 

[22 x 1 — x J\ K 7r x l—x J 

with the momentum fraction, x — fc + /fc+, and < x < 1. 

Introducing in Eq. (6) the relative transverse qq momentum, as defined in Ref. [36], 

K L = (l-a;)k 1 -a;(k T -k)x ) (7) 

we have 

. 2 x 1 f dxdrK± 1 

(q ' = 2(2^ J xjl - x) ((K - xq)l + M 2 ){ml - M 2 ) ' () 

where is the pion mass. The free-mass operator for qq pair is written in terms of the momentum 
fraction, x, and the relative perpendicular momentum as: 

The matrix element of the neutral pion decay driven by the Abelian anomaly is given by one 
unique CPT-invariant Lorentz structure. When one of the photons is off-shell the same matrix element 
describes the transition amplitude 7*7 — > ir°, 

(n (K)\j;°\ 7 (k 2 )) = e 2 e^ a pe;q a k^F 77r o(q 2 ) , (10) 

where e!^ is the polarization of the real photon. Using Eqs. (2), (3) and (8), the transition form factor 
is given in our model by, 



a N c M 2 f dxfKx 1 

> 6n3 u J x{ i _ x) ((K - xq )l + M 2 ){ml - M 2 ) ■ 

In the soft pion limit, following Ref. [41], the form factor becomes 



F ^' {0) = 4^U (12) 
and the neutral pion decay width is given by [41], 

r° = ^F* o (0), (13) 



where a is the fine structure constant. A charge radius, (r 2 ) , can be defined via, 



2 dF^ Q {g 2 ) 

T-a — — 



dg 2 



<j 2 =0 



which in the soft-pion limit behaves as [42], 



v/2M 



(14) 



(15) 



To complete this section on form factors, we recall the long-standing result on the light front for 
the elastic pion form factor [7]. Its expression is likewise obtained following the steps of Eqs. (4) to (8) 
with the definition of the elastic form factor (q = p' — p): 



T,(p,p<) = (n+(p')\j; |7r+(p)) = eF„(q 2 )(p + p'), 
The clastic form factor is thus given by 



fAq 2 ) = 



2N r M 2 f dxd 2 K ± M o 



(l-a)qj.-Kj. 
K^+A/ 2 



(2tt)3 fl J X (l - x) {ml - M' 2 ){jnl - M 2 ) ' 



(16) 



(17) 



where the free mass operators, Mq(Kj_) and Mq 2 (K' 2 ), were defined in Eq. (9) in terms of x and 
the relative perpendicular momenta, = (1 — x)(p — k)j_ — xk± and K' ± = + (1 — x)q±- The 
expression for the rms charge radius, in the soft-pion limit, 



V3 1 

2^U 



(18) 



is know as the Tarrach relation [43] . 



2.3 Wave function models 

We can identify an asymptotic pion wave function in Eqs. (11) and (17) and proceed to model one 
which reproduces soft QCD at low momentum transfer and hard perturbative effects for large q 2 . We 
follow Ref. [7] where the following replacement is in order, 

and the wave- function <P(K 2 ) is normalized to one: 



j d 3 K$l(K 2 



(20) 



This is the normalization condition within the framework of Hamiltonian light-front dynamics [37; 44] . 
The functional dependence of K 2 on Ki and x originates in the free mass operator, M§, 

K 2 (K ± ;x) = ^-M 2 . (21) 

From Eq. (11) and Eqs. (19) to (21), we obtain the integral expression for the neutral pion transition 
form factor on the light front, 

_ VKAI f dxd 2 K ± <P„(K 2 ) 
^' {q ) ~ 67rl J x(l - x)^AT (K - xq)l + M 2 ' [ZZ) 

In the limit Q 2 = —q 2 — > oo, pQCD predicts the form factor decreases asymptotically as ~ q~ 2 with 
a limiting value, Q 2 F 1 ^o(Q 2 ) = 2/ ff [4; 5]. 
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Table 1 The model's length scale parameter, r nr [Eqs. (26) and (27)], as a function of the constituent quark 
masses and for f„ = 92.4 MeV given by Eq. (23). The corresponding charge radii are listed next to r nr for both 
the neutral and charged pion. 



Model 


m u , d [GeV] 


r nr [fm] 


< rl [fm] 


<rlo> 1/2 [fm] 


Gaussian 


0.220 


0.345 


0.637 


0.683 




0.330 


0.472 


0.655 


0.552 


Hydrogen 


0.220 


0.593 


0.795 


0.782 




0.330 


0.708 


0.807 


0.582 


Experiment [28] 






0.672 ± 0.008 





Table 2 The models' length scale parameter, r nr [Eqs. (26) and (27)], and corresponding pion charge radii for 
increasing values of the pion decay constant f^. 



Model 


/. [MeV] 


"Vd [GeV] 


r nr [fm] 


< rv >V» [f m ] 


< rvo > 1/2 [fin] 


Gaussian 


92.4 


0.220 


0.345 


0.637 


0.683 




97.0 


0.220 


0.303 


0.589 


0.657 




110.0 


0.220 


0.172 


0.406 


0.564 


Hydrogen 


92.4 


0.220 


0.593 


0.795 


0.782 




97.0 


0.220 


0.543 


0.750 


0.767 




110.0 


0.220 


0.410 


0.626 


0.720 


Experiment [28] 








0.672 ± 0.008 





For the purpose of completeness, we also quote the result for the neutral pion decay, k\ = fe| 
The amplitude for this process can be formulated using the same reasoning that led to Eq. (22), 



0. 



2 y/N c M 

3 7T 3 / 2 



dxd 2 K ± 



c(l - x)y/ME ((K - 2xk n 



M 2 ) 



$.{K 2 ) 



where |fc 7 | = m w /2 and x < 1/2. In the soft-pion limit and for a constant pion vertex, <P n 
Mo 2 )- 1 , 



F J7r o (Q 



0) 



(23) 

(-< + 
(24) 



In practice, the calculation of the neutral pion width using either i^c^O) in Eq. (22) or i* 1 7r o_ >77 
in Eq. (23) yield the same numerical result within the light-front approach. This follows from the 
observation that the pion mass is much smaller than the inverse of the characteristic length or mass 
scale of the wave function. 

As mentioned earlier, both form factors are calculated with two different qq bound-state wave 
functions: the Gaussian and the hydrogen-atom model. Both models depend on two parameters, the 
constituent quark mass, M and the length (or mass) scale, r nr , of the wave function. In order to test 
the model's mass dependence, all calculations are done with two constituent-quark masses for which 
we choose M = m U: d = 220 MeV and 330 MeV. This length (or mass) scale is set by fitting the pion 
decay constant, f n '= 92.4 MeV [9], 



<N r M 



dxd 2 k\ 



471- 3 / 2 J x(l - x) v / Mc 
where the two models for the pion wave function are explicitly given by, 

4 



^(K 2 (K ± ;x)) = C^exp 
<P^(K 2 (K ± ;x)) = C v 



-r 2 K 2 
3 111 



1 



(25) 

(26) 
(27) 
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Fig. 2 The space-like 7*7 — > it transition form factor, F 17T (Q 2 = —q 2 ), for two different constituent quark 
masses m q = m u ,d — M and f w — 92.4 MeV; see Eqs. (26) and (27) for the two model wave functions and 
Table 1 for the corresponding length scale parameter, r nr , and charge radii. Note that the ordinate is scaled 
logarithmically. Data are from Refs. [1; 3; 6] . 



where C,r is the overall normalization. Once r nr is fixed, the charge radius (14), the transition form 
factor (22), and the elastic form factor (17) are all calculated with the same given pion wave function 
on the light cone. In Section 3, we explore the model sensitivity to the constituent quark mass with 
regard to the pion's static observables and elastic and transition form factors. In Section 4 we study 
the impact of modifying the pion's length scale r nr on the functional behavior of F 17r o{Q 2 ) and the 
consequences for F 7l (Q 2 ) and the rms charge radius. 



3 Numerical results and variation of the constituent quark mass 

Numerical results for the transition form factor, F J7r o(Q 2 ), and its weighted description, Q 2 F 17T a(Q 2 ), 
are presented in Figs. 2 and 3 for the choice M = m Ut d = 220 MeV and 330 MeV of constituent quark 
masses. The length scale, r nr , is adjusted in both cases so as to reproduce the pion decay constant in 
Eq. (25). As seen from Table 1, the calculated rms charge radius is in agreement (but too small) with 
the experimental value for the Gaussian case but about 20% larger when using the hydrogen model. In 
the asymptotic limit, pQCD predicts an upper bound lin^^,^ Q 2 F 1 ^o(Q 2 ) = 2f n = 0.185 GeV; Fig. 3 
gives the impression that in all but one case the curves are not bounded by this prediction. However, we 
checked numerically that for large enough values of Q 2 all curves converge to limQ2^ oc Q 2 F lv o{Q 2 ) = 
0.2 GeV ~ 2/^. We ascribe the difference of about 8% between this limit and that of pQCD to our 
modeling of the transversal component of the pion wave function. 

Whilst the hydrogen bound-state model for m u _d = 220 MeV yields a form factor, Q 2 F 1 ^o (Q 2 ), 
that comes closest to the BaBar data [1], it does not provide the logarithmic form these new data 
points seem to imply for Q 2 > 15 GeV 2 . Moreover, for Q 2 < 10 GeV 2 , the hydrogen bound-state 
model reproduces very poorly older CLEO data. If one does not take into account BaBar data points 
at large four-momentum-squared transfer, the Gaussian wave function with m u ^ = 220 MeV provides 
the best model description of all data including Belle's measurements and coincides with the pQCD 
limit for large Q 2 . 
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Babar (2009) 
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Q 2 [GeV/cf 



30 



40 



Fig. 3 The momentum-squared weighted transition form factor, Q 2 F- TK {Q 2 ), for two different constituent 
quark masses m q — m Ui d = M and = 92.4 MeV. The dotted horizontal line is the pQCD prediction [5]; see 
Eqs. (26) and (27) for the two model wave functions and Table 1 for the corresponding length scale parameter, 
r nr , and charge radii. Data are from Refs. [1; 3; 6]. N.B. For Q 2 — >• oo, Q 2 F J7r (Q 2 ) reaches the asymptotic limit 
2/tt — 0.185 GeV for both models and masses. 



The elastic form factor for two constituent quark masses and both models is depicted in Figs. 4 
and 5. Here, both the hydrogen and Gaussian models with a constituent mass of 220 MeV are in 
good agreement with data whereas for m Ut d = 330 MeV the form factor appears to be too soft. The 
hydrogen wave function for m Ut d = 220 MeV is particularly successful in reproducing the 1/Q 2 tail 
of Fn(Q 2 ). However, as just noted, the resulting charge radius (r 2 ) 1 / 2 = 0.795 fm is 20% larger than 
the experimental value in this case. On the other hand, this may be a desirable aspect of the model 
since it is crudely given by a quark- antiquark core with a constant dressed quark mass. It is known, 
however, that dressing the quark-antiquark scattering kernel with pion loops introduces an attractive 
force which results in a decrease of the bound state's mass and charge radius (see, e.g., the discussion 
about pion loop effects in Ref. [45]). 

We have tested the present light-front model for various constituent masses from 180 MeV to 
330 MeV and find the mass range 220-250 MeV to be the most consistent with a realistic description 
of the observablcs discussed herein (for a recent and more detailed discussion, see Ref. [15]). While we 
acknowledge the mass dependence in our calculations, this is hardly limited to the present light-front 
model [46]. Indeed, it is a recurrent and unavoidable feature of all constituent quark models which do 
not account for dynamical chiral symmetry breaking. The use of a constituent quark mass is justified 
and yields robust results in calculations which involve heavy-quark propagation [47; 48; 49], however 
this is not true for observablcs involving light quarks [30; 31]. 

Nonetheless, this study aims at a consistent treatment of all form factors, the charge radii and 
decay constant with the same single parameter set, r nr and M = m u ,d- We thus employ in the following 
rn>u d = 220 MeV which reproduces well F^+(Q 2 ), (r„) and /„. and yields a satisfying description of 
FjUQ 2 ) below Q 2 ~ 15 GeV 2 . 
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Q 2 [GeV/cf 

Fig. 4 The space- like elastic form factor, JV(<2 2 ), for two different constituent quark masses m q = m u = 
rrid = M and /•„- = 92.4 MeV; see Eqs. (26) and (27) for the two model wave functions and Table 1 for the 
corresponding length scale parameter, r nr . Data are from Ref. [50; 51; 52; 53]. 



4 BaBar data, the pion's length scale and the impact on the charge radius 

Having fixed mass and length scale parameters of the form factors derived in Section 2, we now allow 
for variation of the latter scale, r nr , while keeping m Ut d = 220 MeV fixed. Table 2 illustrates that when 
r nr decreases so does the rms charge radius, (r 2 ) 1 / 2 , whereas the weak decay constant, / w , increases 
in agreement with relation (18). This is observed independently of the wave function model employed. 
Bearing this on mind, we modify r nr , in other words the pion width in momentum space, and observe 
the impact on the form factors discussed in Section 3. 

The results for the transition and elastic form factors weighed by Q 2 are shown for two different 
charge radii and corresponding decay constants (see Table 2) in Figs. 6 and 7 and for comparison, 
we plot the form factors for both models. The comparison within Fig. 6 and between Figs. 3 and 6 
highlights that the model calculation of Q 2 F lv o(Q 2 ) coincide increasingly better with BaBar data for 
smaller charge radii. In particular, the hydrogen type bound-state function for which = 97 MeV, viz. 
about 5% larger than in Table 1, yields a transition form factor in reasonable agreement with BaBar 
data for Q 2 > 15 [GeV/c] 2 . Nonetheless, below Q 2 ~ 15 [GeV/c] 2 , this model produces a transition 
form factor much harder than experimentally evidenced. For the larger decay constant, /„■ = 110 MeV, 
this is even more so the case. At any rate, the wave function model is not appropriate since it neither 
reproduces the extent CLEO data nor the BaBar data which hints at a logarithmically increasing form 
factor at large Q 2 . 

The best result depicted in Fig. 6 is given by the dashed (blue) line for the Gaussian model and 
with f n = 110 MeV from which we deduce (r 2 ) 1 / 2 = 0.406 [fm], i.e. a pion charge radius 40% smaller 
than its experimental value [28] . Decreasing even more r nr will lead to a harder form factor F J7T o (Q 2 ) 
and improve the agreement with BaBar data, however the corresponding weak decay constant and 
charge radius turn out to take unrealistic values. Moreover, within the light-cone framework, a wave- 
function model congruent with the apparent logarithmic increase of BaBar's F J7T o(Q 2 ) data at large 
Q 2 , e.g. the dashed (blue) line in Fig. 6, is shown in Fig. 7 to fail in reasonably describing the elastic 
form factor, F^IQ 2 ). The hydrogen-type wave functions yield reasonable rms charge radii and Ft T {Q 2 ) 
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Q 2 [GeV/cf 

Fig. 5 The momentum-squared weighted elastic form factor, Q 2 F W (Q 2 ), for two different constituent quark 
masses m q = m u — md = M and / T = 92.4 MeV; see Eqs. (26) and (27) for the two model wave functions and 
Table 1 for the corresponding length scale parameter, r nr . Data are from Ref. [50; 51; 52; 53]. 



form factors yet, as noted before, the transition form factors are poorly described in most of the Q 2 
domain. Unsurprisingly, here too smaller charge radii are the results of harder form factors which are 
in disagreement with the current set of data [50; 51; 52; 53] on the elastic pion form factor. 



5 Conclusive Remarks 

A scries of recent articles [16; 17; 18; 19; 20; 21; 22; 23; 24; 31; 54; 55; 56; 57; 58; 59; 60] have dealt with 
the discrepancy between BaBar data [1] on 7*7 — ► ir° and a long-standing prediction of pQCD [4; 5]. 
Attempts to reproduce the BaBar data can be generally identified within two broad classes: 

• alteration of the asymptotic form of the pion's distribution amplitude or wave function; 

• dressing of the jqq vertex with phenomcnological interactions, such as vector dominance. 

As has been discussed extensively in Ref. [31; 55; 56; 57; 58; 59] the modifications of the pion 
distribution amplitude, as proposed for instance in Refs. [16; 17; 19; 20; 21; 22], deviate drastically 
from its QCD asymptotic form. These altered distributions, 4>(x), are constant or at least non- vanishing 
for x = 0, 1 and characterize an essentially point-like pion. One comes to a similar conclusion within 
the framework of a constituent quark model on the light cone or, for the matter, with any other 
quark model. A faithful reproduction of the BaBar data at large-Q 2 values leads to a non-vanishing 
asymptotic pion wave function in momentum space 1 , which we showed to correspond to an excessively 
small pion charge radius. More generally, if the bound-state wave function does not asymptotically fall 
off as (at least) k~ 2 , the elastic form factor, (Q 2 ), becomes harder and may even turn out to take a 
constant value at large Q 2 . This contradicts the bulk of existing data on the elastic pion form factor. 
We conclude that it is impossible, within a given consistent approach, to describe /„•, (r 2 .) 1 / 2 , F 7T (Q 2 ) 
and F 17T (Q 2 ) equally well. One may alter the pion bound-state function to account for the large-Q 2 
BaBar data; this, however conflicts with other data and QCD-based studies that produce soft pions, 

1 at least for practical purposes, as we use Gaussian and hydrogen wave function models which eventually 
do vanish at large enough k 2 even for extreme values of r nr . 
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Fig. 6 The momentum-squared weighted transition form factor, Q 2 F 11T {Q 2 ), for m q = m u = rrid = M 
220 MeV and two values of f n > 92.4 MeV. See Table 2 for the corresponding length scale parameter, r„ 
Data are from Refs. [1; 3; 6]. 



viz. pion distribution amplitudes that vanish as ~ (1 — x) 2 for x ~ 1. We also note that our pion 
transition form factor, the solid line in Fig. 6, is in very good agreement with predictions by Agaev et 
al. [54], Bakulev et al. [57] and Brodsky et al. [58; 59]. 

With respect to the photon vertex modification [18; 23; 24], whilst resonances dominate the electro- 
magnetic form factor of the process 7* — > 7r + 7r - in the time-like region, this is not necessarily true in 
the deep space-like domain. It is reasonable to assume the imaginary part of the form factor vanishes 
at large time-like momentum squared, s = q 2 — > 00, and has a unitary cut on the real axis starting at 
the threshold s = 4m \. In this case, the real part is obtained from an unsubtracted dispersion relation, 

f , w= ir^, (28) 

7T ,/ 4m 2 S - S' 

The imaginary part receives contributions from intermediate resonances that dress the quark-photon 
vertex. Most of these intermediate states, 7* — > X — > 7r + 7r~, lie below s = 10 GeV 2 , the most prominent 
of which being the p°(770). However, the integration in Eq. (28) is over infinite values of s' and thus, in 
averaging out over all possible intermediate states, these resonances are unlikely to have a big impact 
on the real part of F 7r (s). 

Indeed, analyticity and unitarity suggest that higher resonances have little impact on the space- 
like elastic form factor even at high momentum transfers, which was shown with a phcnomcnological 
light-front constituent quark model that incorporates photon dressing [12; 13]. This model describes the 
space-like region up to —10 [GeV/c] 2 , while in time-like region results are reasonable up to 10 [GeV/c] 2 . 
The approach was developed to calculate the elastic pion form factor in the space- and time-like 
regions starting from the covariant Mandclstam formula. A vector meson dominance model, built 
microscopically from the resolvent of a light front mass-squared operator in the spin-1 channel, was 
used for the quark-photon vertex dressing. This resolvent is obtained from the vector meson resonance 
wave functions. The wave functions of both the pion and vector meson resonances are eigenstates of 
the relativistic constituent quark squared-mass operator [61; 62], which accounts for both confinement 
through a harmonic oscillator potential and ir ~ p splitting by means of a Dirac-delta interaction in the 
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Fig. 7 The momentum-squared weighted elastic form factor, Q 2 F n (Q 2 ), for m q = m u = = 
and two values of > 92.4 MeV. See Table 2 for the corresponding length scale parameter, r B 
Ref. [50; 51; 52; 53]. 



M = 220 MeV 
.Data are from 



pscudoscalar channel. The coupling of the vector resonances to the photon and the decay vertex of the 
vector meson to the quark- antiquark pair are calculated within the model. The vector decay constants 
and the overlap between the vector-meson resonances and pion light-front wave functions decrease fast 
with the radial excitation. This is because the increasing number of nodes in the excited-state wave 
functions of vector mesons naturally damps both the decay constant as well as the overlap with the 
pion, strongly suppressing the contribution of the higher resonances to the space-like pion form factor. 
This behavior can be appreciated in Fig. 9 of Ref. [13] where it is seen that the higher resonances have 
no influence on the functional behavior of the elastic form factor. 

If an analogous model is applied to the 77* — > tt° transition form factor, one expects the higher 
resonances to be suppressed as elucidated by the general discussion based on the dispersion relations 
in Eq. (28). There is no obvious sign that any amount of dressing the quark-photon vertex with higher 
resonances can significantly alter the asymptotic q~ 2 decrease of the form factor, whether in the region 
Q 2 = 20 — 40 [GeV/c] 2 or below. Of course, this issue merits further investigation, which we postpone 
to future scrutiny. 

The results of our study suggest the large-Q 2 BaBar data are inconsistent with other experimental 
data on form factors, in particular the Belle data, and static properties of the pion. These are, in turn, 
in agreement with long-standing numerical results of pQCD and soft QCD. Moreover, it is intriguing 
that the 7* — > 777 and 7* — > 7/7 transition form factors also measured by the BaBar Collaboration at 
Q 2 = 112 [GeV/c] 2 [63] are in full agreement with the CLEO results [3] and theoretical results, see for 
instance the green band in Fig. 2 of Ref. [57]. 
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